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ABSTRACT 

In a companion paper (Arnouts et al. 2004) we presented new measurements of the galaxy luminosity 
function at 1500A out to z~l using GALEX- VVDS observations (1039 galaxies with NUV<24.5 and 
z>0.2) and at higher z using existing data sets. In this paper we use the same sample to study evolution 
of the FUV luminosity density pisoo- We detect evolution consistent with a (l+z) 2 - 5±0 - 7 rise to z^l 
and (l+z) 05±0 4 for z>l. The luminosity density from the most UV-luminous galaxies (UVLG) 
is undergoing dramatic evolution (x30) between 0<z<l. UVLGs are responsible for a significant 
fraction (>25%) of the total FUV luminosity density at z~l. We measure dust attenuation and star 
formation rates of our sample galaxies and determine the star formation rate density (p*) as a function 
of redshift, both uncorrected and corrected for dust. We find good agreement with other measures of 
p\ in the rest ultraviolet and Ha given the still significant uncertainties in the attenuation correction. 
Subject headings: galaxies: luminosity function — galaxies: evolution — ultraviolet: galaxies - 
observations: cosmology 



1. INTRODUCTION 
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The rest-frame far-ultraviolet (FUV; 1500A) luminos- 
ity has been used to determine the star formation rate 
(SFR) of stellar populations over the complete range of 
rcdshifts for which galaxies have been observed. The 
utility and limitations of the integrated measures — the 
FUV luminosity function (4>fuv) an d luminosity den- 
sity (/?i5oo) — and their relation to the star formation his- 
tory of the universe has been extensively discussed and 
reviewed (e.g. Madau, Pozzetti, & Dickinson 1998, Hop- 
kins 2004). A principal goal of the Galaxy Evolution Ex- 
plorer (GALEX) mission (Martin et al. 2004) is to per- 
form deep wide-angle surveys to obtain an accurate mea- 
surement of the evolution of the FUV luminosity density 
over the range < z < 1 and beyond. In this letter we 
present results from a small pilot study performed in the 
2h VIRMOS-VLT Deep Survey field using measurements 
from 1039 galaxies. 

GALEX data will allow us to determine how the rest- 
UV can best be used to study the detailed properties of 
galaxies (e.g. dust, metallicity, star formation history). 
Since this is work in progress, here we will instead use 
existing methods to determine the intrinsic luminosity 
of galaxies in the FUV (Meurer, Heckman, & Calzetti 
1999, hereafter MHC99, for dust corrections) and the 
SFR that this luminosity implies (Kennicutt 1998 for 
SFR conversion) . This simple analysis yields some quick 
answers — we will discuss how this work will be expanded 
and developed in the near future. 

Throughout this paper we adopt the flat-lambda cos- 
mology (h M = 0.3, Q A = 0.7) with H = 70 km s" 1 
Mpc" 1 . 

2. DATA 
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TABLE 1 
FUV 1500A Luminosity Density 
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a units: erg s 1 Hz 1 Mpc 3 , flat-lambda cosmology with Ho = 
70 km s" 1 Mpc" 1 (n M = 0.3,fi A = 0.7) 
h L min = 0.2L*(z = 3) d ;M min = -19.32 
c data from Wyder et al. (2004) 
d data from Steidel et al. (1999) 
e fixcd faint end slope a = —1.6 
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Fig. 1. — FUV luminosity density vs. redshift. Filled circles 
indicate LF fit to full sample with unconstrained slope a. Filled 
squares denote LF fit with fixed a = —1.6. Purple (GALEX) 
and dark green (HDF) symbols are from this work. Black dot is 
taken from local LF (Wyder et al. 2004) and green square from 
Steidel et al. (1999). Open circles denote pisoo determined using 
~Vmax ■ Errors do not include cosmic variance. Lines indicate (1 + 
z) n evolution. Dotted, solid and dashed lines correspond to n = 
1.5,2.5,3.5 respectively. 



GALEX observations of the VVDS 0226-04 field 
(02h26m00s -04°30'00", J2000) were performed in 
October-November 2004 as part of the GALEX Deep 
Imaging Survey. Further details of these observations, 
the subsequent match to VVDS spectroscopy and pho- 
tometry and the calculation of the LF can be found in 
the companion letter Arnouts et al. (2004) (hereafter Pa- 
per I) and references therein. Paper I also describes the 
derivation of the 4>fuv a t z=2.0 and 2.9 using an HDF 
sample from Arnouts et al. (2002). For comparison we 
also use the local <\>fvv (Wyder et al. 2004) and the z~3 
Lyman Break Galaxy (LBG) (f>FUV (Steidel et al. 1999) 

3. LUMINOSITY DENSITY 

We calculated the FUV luminosity density pi5oo(z) 
from the GALEX-VVDS sample in four redshift bins 
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Fig. 2. — FUV luminosity density of ultraviolet luminous galax- 
ies (UVLGs) vs. redshift and comparison w/ QSO luminosity den- 
sity. Filled circles from P1500 UVLG integrated from 00 down to 
0.2L* :Z — 3 or yipjjv -19.32. Colors same as Figure 1. Vertical 
hatched bars indicate fraction of luminosity emitted by galaxies 
brighter than 0.2L J . jZ=3 . Red dashed line shows QSO FUV LD 
using values from Boyle et al. (2000) and Madau et al. (1999). 
Solid line same as Figure 1. 



(< z >= 0.3,0.5,0.7,1.0), and also determined values 
for the HDF sample (< z >= 2.0,2.9). Results are 
shown in Table 1 and plotted in Figure 1. We chose 
to calculate pi5oo( z ) m several ways. First we summed 
4>(L)LdL using the LF obtained from the V max method. 
Because luminosity bins with no detections do not con- 
tribute we consider this a lower limit on pi5oa(z). We 
also calculated a "total" luminosity density by integrat- 
ing Schechter function fits to the LF using the formula: 

dLLcj)(L) 

with L min = (p = 4> te L^,T{a + 2)). Although this quan- 
tity is strongly dependent on uncertainties in the faint 
end slope (a), it allows direct comparison with other 
measurements of pisoo (z) and the star formation rate 
density, p\. Fits and errors were determined using the 
ALF tool (Ilbert et al., 2004) with error bars based on 
the extreme values of the LD calculated at each point 
on the a-M* lcr error contour. For the z=1.0 bin, our 
best Schechter function fit yielded large errrors for the 
slope (a=-l. 63^0 43)- For this bin we fixed the faint end 
slope at a to -1.6, adopting the value used in high-z stud- 
ies (e.g. Steidel et al. 1999) and consistent within errors 
with our own values at lower and higher z. Total pi5oo(^) 
shows significant ~ (1 + z) 2 ' 5 evolution out to z~l, with 
evidence for a shallow continued rise out to z=3. This 
evolution is discussed further in the next section. Two 
points are worth noting regarding the comparison of LD 
at different redshifts. First, as demonstrated in (Paper 
I) and discussed below, the galaxy population that con- 
tributes most of the LD varies (vs. color, luminosity) 
with redshift. Secondly, while most of the sample is UV- 
selected, the Steidel et al. (1999) LBG galaxies were 
color-selected and the z=3 LD value may be missing 
some fraction of the UV light. The similarity between 
the z=2.9 and z=3 data points suggests that the missing 
fraction is small. 
We explore the contribution to the luminosity density 



Evolution of the Far-Ultraviolet Luminosity Density 



3 



Afuv 

-10 12 




Fig. 3. — Histogram of FUV slope Pglx f° r 888 galaxies with 
measurements in U-band (black) and split into two subsamples: 
2 <0.5 (brown) and z>0.5 (orange). (Red/white) Distribution 
of /3glx for "isolated" GALEX detections with only one optical 
counterpart within 4" radius. 
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from UV luminous galaxies (UVLGs) by measuring the 
luminosity density from galaxies with L > L m i n . To 
facilitate comparison with high-z studies, we set L m i„ = 
0.2L*, Z=3 (M min = -19.32) from Steidel et al. (1999), 
also adopted by Giavalisco et al. (2004) for their work. 21 
These galaxies are observable in all redshift ranges and 
therefore there is no additional uncertainty related to 
extrapolation beyond the faintest observed magnitude. 
Figure 2 highlights the dramatic evolution of Pi5oo,uvlg, 
increasing by x30 to z^l or (1+z) 5 . Furthermore we find 
that UVLGs are a major contributor to P1500 at z~l with 
a fractional contribution, pi5oo,uv lg / P1500 of 25%. We 
plot for comparison pi5oo,QSO using the functional form 
of the QSO LD evolution (in the B-band) from Boyle 
et al. (2000) and the QSO SED from Madau, Haardt, 
& Rees (1999) which has a shallower evolutionary slope 
(3.5) vs. UVLGs (5) for z<l. 

4. STAR FORMATION RATE DENSITY 

To determine intrinsic ultraviolet luminosities for the 
GALEX- VVDS sample, we apply the MHC99 dust at- 
tenuation formula: 

A FUV = 4.43 + 1.99(/3) = 4.49 + 1.97(p GLX ) 

where we use the definition of (3glx, the FUV slope 
calculated using the rest-frame GALEX FUV and NUV 
bands, from Kong et al. (2004). We only calculate Pglx 
for the subset of galaxies observed in U-band (888 galax- 
ies). Typical errors are 073 ~ 0.4. Figure 3 shows the 
distribution of the k-corrected (3qlx- The full sample 
has median /3 GLX = -1.64, FWHM{(3) = 1.4 with lit- 
tle variation with redshift. We find good agreement with 
measurements of (3 at low-z (< (3 >= —1.6 for a FUV se- 
lected sample; Treyer et al. 2004) and high-z (Adelberger 
& Steidel 2000). 

Within our own sample we might have expected to 
see an increase of Pglx vs. z since high luminosity 
galaxies — which dominate the high-z bins — are expected 
to show significant attenuation. Several effects could 

21 This luminosity corresponds to 10 101 Lq, ~ 2/3 the luminos- 
ity limit (10 10 - 3 L ) adopted for UVLGs in Heckman et al (2004). 



Fig. 4. — Star formation rate of GALEX- VVDS galaxies vs. 
redshift using Kennicutt (1998) conversion. Top: Uncorrected 
star formation rate for galaxies with varying Apjjy Solid line 
and dashed line correspond to NUV<24.5 (current sample) and 
NUV<26 (GALEX Ultra-deep survey) limits. Stars show values 
for L* from Paper I and Wyder et al. (2004). Dotted line corre- 
sponds to L m i n = 0.2L J . jZ= 3 cut. Bottom: attenuation-corrected 
star formation rates. Symbols same as top. Colored lines corre- 
spond to detection limits for NUV<24.5 at minimum attenuation 
level for each subsample. 



work against this trend. We are detecting galaxies close 
to the NUV band confusion limit (beam/source~10 for 
NUV<25) and source blending could shift UV-optical 
colors and the slope blueward. We performed tests which 
conservatively apportioned NUV flux among all poten- 
tial optical counterparts and set a limit on the offset 
of the median AflcLXMend < 0.35. This is consistent 
with the median /3glx = —1-44 measured for "isolated" 
UV detections with only a single optical counterpart (see 
Figure 3). (However, we can't neglect the possibility 
that some fraction with multiple counterparts are phys- 
ical pairs which could show a different distribution of 
Pglx)- We also note that the MHC99 A FUV -f3 relation 
was determined for starbursting galaxies (the bulk of our 
sample, see Paper I) but might overestimate the correc- 
tion for normal star forming galaxies (Bell 2002; Kong et 
al. 2004) which are found in our lowest redshift bins. For 
a conservative measurement of the average attenuation in 
our whole sample, we use the "isolated" subsample, and 
calculate a mean attenuation factor of x7 (A™g^ = 1.8) 
where we have estimated and applied a bias correction to 
the mean (x0.7) due to non-negligible a p. We also adopt 
a 'minimum attenuation' A™jyy = l which may be more 
representative of a full UV-selected population (Buat et 
al. 2004). 

The SFR was calculated for each galaxy using: 

— 1\ i a w in— 28 r /„ — , „— 1 JJ Z ~ 1> 



= 1.4 x 10- 28 Lfw(erg s" 



SFR{M yr" 

from Kennicutt (1998). In Figure 4 we plot the SFR 
derived for each galaxy using the uncorrected and the 
dust-corrected FUV luminosities. Our sample shows no 
dependence of dust attenuation with SFR uncor and as a 



4 



Schiminovich et al. 



-l.o 




go -2.0 : $ 



0.0 



0.5 1.0 1.5 2.0 2.5 3.0 

Redshift z 



Fig. 5. — Star formation rate density vs. z. Filled circles from 
measurements at 1500A (uncorrected for dust) same as in Fig- 
ure 1. Blue comparison points are rest-frame UV measurements 
uncorrected for dust attenuation. Inverted blue triangle from Sul- 
livan et al. 2000. Dark blue triangles from Lilly (1996). Light 
blue triangles from Wilson et al. (2002) for a = —1.5. Solid line 
rises as (1+z) 2,5 for z<l and then (l+z) 5 for z>l based on chi- 
squared fit to our sample (see inset; lcr and 2a confidence contours 
shown). Shaded region shows range corresponding to max/min 
dust-attenuation. Filled red stars from dust-corrected Ha mea- 
surements (with increasing redshift) from Perez-Gonzalez et al. 
(2003), Gronwall (1999), Tresse & Maddox (1998), Tresse et al. 
(2002). Open red star from SDSS (Ha/emission line) Brinchmann 
et al. 2004. 



consequence we find higher attenuation in galaxies with 
high SFR cor . This paucity of low- attenuation galax- 
ies with high SFR cor has been noted in previous stud- 
ies (e.g. Wang & Heckman 1996, Adelberger & Stei- 
del 2000). Some of the observed effect may also be 
due to the scatter in Apuv discussed above (resulting 
in a tail of high Apuv galaxies) and/or limitations of 
the dust attenuation law. We plot p\(z) (derived from 
P1500 with no dust correction) in Figure 5. Measure- 
ments from this paper were fit using the parametrization 
from Baldry et al. (2002) (p*(z) - (1 + zY™° l ,z < 1 
and p\(z) ~ (1 + z) a "^',z > 1). We find a best-fit 
P evo i = 2.5 ± 0.7, cw = 0.5 ± 0.4. The lcr con- 
straint on the {oievoh Pevol) pair is consistent with in- 
dependent derivations using 2dF (Baldry et al. 2002), 
SDSS (Brinchmann et al. 2003) and other recent studies 
(c.f. Fig. 13 Baldry et al. 2002). 

Several uncorrected (blue) and dust-corrected (red) 
comparison measurements obtained using spectroscopic 
redshifts are shown in Figure 5. Before determining p\ 
we converted P2000 (Sullivan et al. 2000, Lilly et al. 
1996) and p 2 500 (Wilson et al. 2002, a = -1.5 data) 



to /O1500 using p(X) obtained from local P1540 and P2300 
by Wyder et al. 2004 (~ A 9 ). Wilson et al. (2002) 
and Lilly et al. (1996) both show good agreement with 
our measured values despite the difference in evolution- 
ary slope obtained in the two studies ((3 evo i ~ 1.7 ± 1, 
3.3 ± 0.7, respectively) . The local luminosity density 
reported by Sullivan (2000) appears high, as noted in 
Wyder et al. (2004). Finally, we show a likely range 
of dust-corrected SFR densities, applying the average 
A™fjy, A™g-y to the best- fit parametrized p\(z). Us- 
ing the Kennicutt (1998) SFR conversion, we find that 
recent dust-corrected Ha measurements fall within our 
attenuation-corrected range. Although we have implic- 
itly assumed no evolution in the dust correction, we em- 
phasize that for UV flux-limited samples we might expect 
evolution in the average dust-attenuation correction vs. 
redshift and will explore this further in future work. 

The FUV is tracing a predominantly homogeneous 
population (star-forming and starbursting) making inter- 
pretation of integrated measures much more straightfor- 
ward than at longer wavelengths (cf. Wolf et al. 2003). 
We have shown that a significant population of UVLGs 
lies within easy reach (0.6 < z < 1.2). We will com- 
pare these unique star-forming galaxies with their high- 
rcdshift LBG analogs (e.g. Shapley et al. 2003). In the 
near future our sample will expand by x5 in this field 
alone, and by more than x 100 using data from redshift 
surveys across the sky. In some locations we will increase 
our depth to m^s ~ 26 as part of the Ultra-Deep Imag- 
ing Survey and probe down to 0.1L* (see Figure 4) to 
better constrain the faint end of 4>fuv ■ This will be sup- 
plemented by an even larger catalog (>10 6 objects) with 
photometric redshifts. We will soon be able determine 
how SFR evolution depends on environment, morphology 
and spectral type and will examine our results within the 
context of cosmological simulations. A major challenge 
lies in the understanding the role of dust obscuration, 
one which we will explore using recent, more sophisti- 
cated models (e.g. Kong et al. 2004) as the GALEX 
surveys continue. 

GALEX (Galaxy Evolution Explorer) is a NASA Small 
Explorer, launched in April 2003. We gratefully acknowl- 
edge NASA's support for construction, operation, and 
science analysis for the GALEX mission, developed in 
cooperation with the Centre National d'Etudes Spatiales 
of France and the Korean Ministry of Science and Tech- 
nology. 
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